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^ ABSTRACT 

t-H We study the far-infrared (IR) and sub-millimeter properties of a sample of ultraviolet 

(UV) selected galaxies at z ~ 1.5. Using stacking at 250, 350 and 500 fim from Herschel 
Space Observatory SPIRE imaging of the COSMOS field obtained within the HerMES 
key program, we derive the mean IR luminosity as a function of both UV luminosity 
^ and slope of the UV continuum f3. The IR to UV luminosity ratio is roughly constant 

V™j over most of the UV luminosity range we explore. We also find that the IR to UV 

luminosity ratio is correlated with f3. We observe a correlation that underestimates the 
correlation derived from low-redshift starburst galaxies, but is in good agreement with 
the correlation derived from local normal star- forming galaxies. Using these results we 
reconstruct the IR luminosity function of our UV-selected sample. This luminosity 
function recovers the IR luminosity functions measured from IR selected samples at 
the faintest luminosities (Lir ~ 10 11 Lq), but might underestimate them at the bright- 
end (Lir > 5 x 10 11 Lq). For galaxies with 10 11 < Lir/Lq < 10 13 , the IR luminosity 
function of a UV selection recovers (given the differences in IR-based estimates) 52-65 
to 89-112 per cent of the star-formation rate density derived from an IR selection. 
The cosmic star-formation rate density derived from this IR luminosity function is 
61-76 to 100-133 per cent of the density derived from IR selections at the same epoch. 
Assuming the latest Herschel results and conservative stacking measurements, we use 
a toy model to fully reproduce the far IR luminosity function from our UV selection 
at z ~ 1.5. This suggests that a sample around 4 magnitudes deeper (i.e. reaching 
u* ~ 30 mag) and a large dispersion of the IR to UV luminosity ratio are required. 

Key words: ultraviolet: galaxies - infrared: galaxies - methods: statistical. 
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1 INTRODUCTION 

Star-formation is one of the main properties used to trace 
galaxy formation and evolution. Our ability to constrain 
the mechanisms that drive galaxy evolution hence depends 
to a large extent on our ability to measure accurate star- 
formation rates for samples of galaxies. A number of star- 
formation tracers are routinely used, from the strengths 
of spectral lines to broad-band measurements ( jKennicutt] 
1998). Broad-band measurements offer the advantage that 



estimates can be derived for a large number of galaxies 
with minimal spectroscopic follow-up, albeit at the expense 
of possible contamination by strong spectral features (e.g. 
Smail et al.|201l" ). One of the main star- formation rate trac- 
ers of this kind comes from the ultraviolet (UV) range of 
the spectrum, where most of the ener gy is emitted by youn g 



stars (ages ~ 10-10 8 Myr; see e.g. Martin et al. 2005a) 



The UV has been widely used over a large redshift range to 
infer the cosmic star-formation density, from z = 7 to z = 
(e.g. |Schiminovich et al.|2005| |Bouwens et~aLl|2009 ). 

Optical/near infrared (IR) observations probe the rest- 
frame UV for samples of galaxies at high redshift, hence 
it has been the primary choice for constraining the star- 
formation activity of the Universe at early epochs from large 
samples. However, interstellar dust, which is a byproduct 
of star-formation, makes the measurement of star-formation 
activity challenging at these wavelengths. Dust grains scat- 
ter or absorb the light emitted by young stars; hence only a 
fraction of the energy output from star-formation is observ- 
able in the UV. Dust grains re-emit this energy over the full 
IR range 8 — 1000 /xm. One way to estimate the dust attenu- 
ation in the UV and to assess the selection bias inherent to 
the UV is then to study the far IR properties of UV-selected 
galaxies. 

A number of studies have used this approach to charac- 
terise the amount of dust attenuation. In the local Universe, 
the star-formation rate density is roughly equally divided be- 
tween UV and IR contributions (Martin et al.|2005b| |Both- 



well et al.|2011| ). At earlier epochs, however, the fraction of 
star-formation rate density, which is directly measurable us- 
ing the UV continuum, decreases from 44 per cent (z — 0) to 
roughly 15 per cent at z ~ 1 ( Takeuchi, Buat, &; Burgarella| 
2005| |Tresse et aL|2007) , w hile it might increase slightly to 
20 per cent at 2 < z < 3.5 ( |Reddy et"aL|2008|), and even to 
higher values in the early Universe ( Bouwens et al.||2010 ). 

To overcome this drawback and the lack of deep IR 
data, it is common to use empirical recipes to correct UV 
for dust attenuation. The most well-known is the relation 
between the slope of the UV continuum and the ratio of 



the luminosities in the IR and the UV (Meurer, Heckman, 



& Calzetti 1999}. The slope of the UV continuum can be 



derived from rest-frame UV colours, and hence is convenient 
for estimating star formation rates at high redshifts when a 



limited wavelength range is available ( |Schiminovich et al. 
~ Bouwens et al.||2009 ). 



2005 



However, this recipe does encounter several pitfalls: it 
has been derived from local starburst galaxies, and it might 
not be valid for more normal sta r-forming galaxies (|Boissier 



et al.||2007| |Cortese et al.||2006| |Munoz-Mateos et al.| 2009 
Hao et al. |2011| ). Moreover, the relation between dust atten- 
uation and UV slope depends on the selection criteria ( |Buat 
et al. 2005; Se ibert et al.|2005| , and might also be sensitive 



to star-formation history (Kong et al. 2004; Panuzzo et al.| 
2007 ), as well as dust properties ( |Inoue et al.|20Q6| ) and dust 
geometry ( |Calzetti|200T| . 

It is hence of particular importance to follow the evolu- 
tion with redshift of the IR properties of UV-selected galax- 
ies, in order to characterize the biases inherent in such se- 
lection, and also to examine the validity of the empirical 
recipes commonly used to correct for dust attenuation. 

In this context a new era starte d with the availabili ty of 
data from the Hersche Q telescope ( |Pilbratt et al.|2010) . In- 
deed, while Spitzer data uncovered the dusty star- format ion 
history of the Universe up to z ~ 1 ( Le Floc'h et al.|2005 ), at 
higher redshifts large extrapolations are needed to estimate 
IR luminosities from 24 /xm data, which could lead to sys- 
tematic errors ( |Bavouzet et al.|2008a| , as they do not probe 
the peak of the dust emission. Based on Herschel data, |E1-| 
baz et al. (2010) showed, for instance, that using mid-IR 



data at z > 1.5 leads to an overestimation of the total IR 
luminosity. Another important feature of the submillimeter 
wavelength range is that the contribution of Active Galactic 
Nuclei to the galaxy Spectral Energy Distribution (SED) is 
generally outweighed by the star-formation component for 
A > 30 fim ( Hatziminao glou et al.|2010 ). 

In this paper, we focus on a UV-selected sample at 
z 1.5 to study with unprecedented statistics the far IR 
properties of UV-selected galaxies using Herschel data for 
more accurate measurements of IR luminosities. Given the 
confusion-limited nature of these data, we rely on a stack- 
ing analysis to derive the mean IR luminosities for different 



classes of object (e.g. Bethermin et al. 
2012||Viero et al.|2012b| ). 



2012 



Hilton et al. 



The paper is organized as follows. We start by present- 
ing the sample construction and the need for stacking (Sec- 
tion. [5]). In Section [3] we present our methods for stacking 
measurements and corrections for biases. Section [4. II shows 
our stacking results as a function of UV luminosity and Sec- 
tion [4?2]as a function of the slope of the UV continuum. In 
Section [5] we reconstruct the total IR luminosity function 
of our UV-selected sample using the stacking results, and 
examine the implications of these results for the measure- 
ment of the cosmic star-formation density from UV and IR 
selected samples. We discuss these results in Section [6] and 
present our conclusions in Section [7] 

In this paper, we use a standard cosmology with Q m = 
0.3, Ov = 0.7 and H = 70 km s -1 Mpc"\ denote FUV and 
IR luminosities as vL v , and use AB magnitudes. 



2 DATA SAMPLE 



We use optical imaging of the COSMOS field from[Capak 
et al.| 02007) in the ?x*-band (obtained at CFHT; depth: 
26.4 mag at 5cr for a 3arcsec aperture) and Vj-band (from 
Subaru; depth: 26.6 mag, also at 5cr for a 3arcsec aperture). 
We generated catalogues from the u* and Vj images us- 
ing SExtractor (Bertin & Arnouts 1996) in order to select 
galaxies directly in the ix*-band, and obtain accurate total 
fluxes. Comparison with the photometry from |Capak et ah] 

1 Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA. 
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(2007) shows good agreement. Hereafter all quoted magni- 



tudes are corrected for Galactic extinction using dust maps 



from Schlegel, Finkbeiner, & Davis (1998) 



We estimated the incompleteness by injecting 1000 fake 
sources in the u* image with random positions. We injected 
sources only within areas not masked for edges or bright 
stars. We assumed that the objects in our UV-selected sam- 
ple are pure exponential disks. We assumed that the disk 
scale length for galaxies with Lfuv = L* ~ 1O 1O L0 is equal 
to 3 kpc. This value is 50 per cent larger than the one from 
ithi et ah] ( |2012| , determined at 0.3 < z < 2.0 in the 
band (restframe wavelength 4600 A at z = 1), in order to 
account for the fact that disks are larger at UV wavelengths. 
The value we assumed is in agreement with predictions from 
|Boissier &; Prantzos| (|2001|, as well as with the U V restframe 
measurements from | Ferguson et al. (2004), taking into ac- 
count that the mean luminosity of their sample is around 
5 x L FUV (z = 1.5). We further assumed that the disk scale 
length varies with luminosity as L 1 / 3 , as observed in the 
local Universe ( de Jong &; Lacey||2000 ). We generated fake 
objects following the observed joint distributions of magni- 
tude, redshift, and minor to major axis ratio; we also allowed 
random position angles. Then for each fake object we could 
infer its UV luminosity and its disk scale length. We in- 
jected these objects in the u* image after convolving them 
with the ?x*-band PSF (FWHM = 0.9 arcsec, |Capak et aL 



2007). We performed source extraction with SExtractor on 



this new image using the same parameters as the ones used 
to generate the it* -band catalogue, but using this time the 
ASSOC mode to cross-match directly the detections with our 
input fake sources list, with a 0.5 arcsec search radius. We 
performed this process 1500 times to have sufficient statis- 
tics to quantify the effects of incompleteness on stacking (see 
Sect. 3.1 ). In order to test the impact of real sources on the 
flux estimation of the fake sources, we cross-matched the 
fake sources we recovered with the full u*-band catalog. We 
then rejected fake sources whose magnitudes are perturbed 
by closeby real sources, based on the angular separation and 
the difference between the input and the recovered magni- 
tude. 

We estimate that the completeness is 76 ± 1 per cent at 
the limiting magnitude u* = 26 mag (Fig.[T]). 

We cross-matched our sample with an updated version 



of the photometric redshift catalogue of|Ilbert et al. (2009 



v.2.0). This version differs from the original |Ilbert et al. 



2009| catalogue by the inclusion of additional near-IR pho- 



tometry (YJHK, McCracken et al.||2012| which improves 
the accuracy of the photometric redshifts, in particular in 
the redshift range we are interested in. 98.5 per cent of the 
sources from our sample have a counterpart in the Ilb ert et| 



al. 



12009J) catalogue within a 1 arcsec search radius. 
We built a rest-frame UV-selected catalogue considering 
objects with 1.2 < z p hot < 1.7 and u* < 26 mag. The aver- 
age photometric redshift accuracy at la in (1 + z) is 0.04 for 
this sample, and the mean redshift is (z p hot) = 1.43. At this 
redshift, the effective wavelength of the t£*-band filter corre- 
sponds to 1609 A, which is in the far-UV (FUV) rest-frame. 
Our final catalogue contains 42,184 objects over 1.68 deg 2 , 
after masking of edges and areas around bright stars. 

We use Herschel-SPIRE (IGriffin et al.|[20Tol ISwin 



yard et al.||2010) imaging at 250 (FWHM 
350 (FWHM = 25.15 arcsec 



18.15 arcsec) 
and 500 /mi (FWHM = 




Figure 1. Completeness in the u* band for the UV-selected cata- 
logue. Error bars are the errors on the mean. The vertical dashed 
line shows the magnitude limit we adopt. 



36.3 arcsec) of the COSMOS field obtained as part of 
the Herschel Multi-Tiered Extragalactic Surve 30(HerMES, 
Oliver et al.|2012| ) programme. We use here the images pro- 
duced by th e SMAP pipeline ( |Levenson et al.|[20To{ |Viero 
et al.|2012a ). The effective area (after removing masked re- 
gions) of the overlap between t£*-band and Herschel-SPIRE 
images is 1.5 deg 2 ; 38,074 galaxies from our UV-selected cat- 
alogue are within this area. 

We performed cross-matching between our UV-selected 
catalogue and the SCAT (Smith et a"T1|2Q12| ) blind detec- 
tions at 250, 350 and 500 /im, using a 5 arcsec search ra- 



dius. Considering as 'detected' only Herschel-SPIRE sources 
with fluxes larger than the confusion limit (at 5cr, 24.0 mJy 
at 250 jum, 27.5 mJy at 350 /xm, and 30.5 mJy at 500 /mi, 
Ngu yen et al.|2010 ), we find that less than 1 per cent of the 
UV sources are detected at 250, 350, and 500 /mi. This re- 
sult implies that we need to use a stacking analysis in order 
to study in a statistical way the IR properties of the UV- 
selected galaxies in our sample. In the following, we include 
all UV-selected galaxies within the HerMES footprint in the 
stacking analysis, whether they are detected at the Herschel- 
SPIRE wavelengths or not. Excluding UV-selected galaxies 
detected at SPIRE wavelengths from the stacking input lists 
does not significantly impact our results. 

We further use this sample of UV-selected galaxies de- 
tected at SPIRE wavelengths to determine the dispersion in 
Lir/Lfuv as a function of Lfuv- This sample has a mean 
infrared luminosity (Lib) = 1O 12 L , is slightly brighter in 
UV than our full sample ((Lfuv) — 2.6 x 1O 1O L0, compared 
to 1O 1O L0 for the full sample), and has a mean IR to UV 
luminosity ratio of (Lir/Lfuv) = 66. 



2 http://hermes.sussex.ac.uk 
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3 STACKING MEASUREMENTS 



We u se the IAS library ( |Bavou zet 2008b ; Bethermin et al. 
|2010a p]to perform the stacking. We use the calibrated 250, 
350 and 500 /im images, and we do not attempt to clean the 
image of any detected sources at SPIRE wavelengths. 

For a given stacking measurement, we generate both a 
postage stamp image and a radial profile, using the mean of 
the individual images that are included in the stacking. We 
derive errors on the profiles by bootstrap resampling. 

The background for the stack images is usually consid- 
ered constant and determined by the average value of the 
stack profile at large distances from the centre of the stack 
image. However, a number of effects can yield a non homo- 
geneous background; we correct our stacking measures for 
two effects which have an impact on the background: the 
incompleteness of the input catalogue, and the clustering of 
the input galaxies. These corrections are based on the work 



of |Bavouzet| ( |2008b| and |Bethermin et al.| ( |2010b ). 



3.1 Correcting for stacking bias 

While limited, the incompleteness at the faint-end of the 
input catalogue can have an impact on the stacking mea- 
surements, for instance if it is partly related to the local 
geometry during the detection process. In particular, the 
detection efficiency for faint objects is lower in dense areas 
of the image. As we will see, this effect can also be described 
by a clustering term under the form of a cross correlation 
between sources of different u* band fluxes. 

If we stack a population of objects which is randomly 
distributed in the sky, we will get a flat background. On the 
other hand, if there is a bias introduced by the detection 
process, we will miss the contribution of the objects which 
are not recovered by the source extraction. The basis of this 
correction is to measure the actual background of the stack- 
ing for a given class of input sources. 

We assume here that the stacking bias effects are re- 
lated to the UV luminosity of the objects. In other words, 
we consider the actual Lpuv distribution for each class of 
galaxy to correct for this bias. To quantify the impact of the 
stacking bias, we use the fake sources created to estimate 
the completeness of our catalogue (see Sect. [2|. We stack 
the fake sources recovered by the detection process with the 
same UV luminosity distribution as the class of galaxy we 
are considering. We show in Fig. [2] the radial profiles of the 
stacking as a function of Lfuv- These profiles are based on 
a sample of fake sources with roughly 25 times the actual 
number of galaxies in our sample. The error bars on these 
profiles are obtained through bootstrap resampling. For faint 
objects (Lfuv < 10 10 Lq), at smaller scales, the profiles are 
actually lower than zero. Note that the amplitude of this 
effect increases for fainter objects. 

This result is related to the well-known effect that the 
detection efficiency is lower in dense areas, and in partic- 
ular for faint objects which are close to brighter ones. If 
all fake sources were to be recovered by the UV source ex- 
traction, the profile of their stacking at the Herschel- SPIRE 



3 |http: //www. ias .u-psud. fr/irgalaxies/f iles/ias_ 
stacking_ lib . tgz| 



wavelength would be zero, as their input distribution is in- 
dependent of that of the real sources. However, the faint 
sources which are close to bright real sources in the u* image 
are not completely recovered. This effect is more important 
for fainter sources, as well as for smaller distances with re- 
spect to bright sources. Hence, the stacking of faint objects 
is missing the contribution to the background of the ones 
which happen to be closer to UV-bright sources, and then 
the background at small scales is lower. 

In order to correct our stacking measurement for this 
effect, we subtract these profiles, without smoothing, from 
the profiles of the stacked images. We add in quadrature the 
errors of these profiles to the errors of the stacked image 
profiles. 

We show as squares in Fig. [3] the ratio of the flux den- 
sities corrected from the stacking bias to the flux densi- 
ties measured with PSF-fitting, assuming a constant back- 
ground, as a function of Lpuv- The amplitude of the cor- 
rection is maximal at faint luminosities (around 2 times 
the flux from direct PSF-fitting at Lfuv ~ 6 x 10 9 Lq) 
then decreases with luminosity, to be negligible for Lfuv > 
2xlO lo L . 

3.2 Correcting for clustering of the input 
catalogue 

We use the formalism developped by |Bavouzet| ( 2008b ) and 
|Bethermin et alT| ( |2010b| in order to take into account the 
impact on the measured flux from the clustering of the pop- 
ulation under study. If the input population is uniformly 
distributed on the sky, the measured stacking is just the av- 
erage flux of the population at the stacked wavelength. In 
practice, galaxies are clustered, so there is an excess of prob- 
ability to find another galaxy of the sample within the beam, 
compared to the value derived from a randomly distributed 
population. This yields an over estimation of the flux. The 
probability is proportional to the angular correlation func- 
tion of the class of galaxies under study. The two dimen- 
sional profile of the resulting stacking can then be written 



1(0, = S x PSF(0, 0) + a x [w(0, 0) * PSF(0, 0)] . (1) 

Here S is the average flux, PSF(0, cj>) is the Point Spread 
Function at the stacked wavelength, w(0,c/)) is the angular 
autocorrelation function of the input population; the symbol 
* denotes a convolution, and a is a parameter relating to the 
density of the input population and their contribution to the 
cosmic background. The effect of the clustering is to add to 
the profile a component which is broader than the pure PSF 
(see e.g. Beth ermin et al.|2012 their fig. 3). 

To correct our stacking measures for clustering, we ad- 
just the radial profile of each stacked image over 120arcsec 
following eq. [I] with S and a as free parameters, and using 
the autocorrelation function w(6) of the sample. We obtain 
the best value of S by marginalising the two dimensional 
probability of (5, a) over a, using % 2 statistics. 

In this paper we perform stacking in bins of UV lumi- 
nosity and the slope of the UV continuum (/?). We measure 
w{6) using the method of Szapudi et al. (|2005|, and fit it 



with a power law w{6) = A w 9 , correcting for the integral 



constraint following Roche & Eales (1999). We find that the 
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Figure 2. Profiles of the stacking at 250, 350 and 500 jim of recovered fake sources in bins of Lpuv luminosity. 



correlation function is well modelled with 5 — 0.64 ± 0.07. 
Note that we consider a as a free parameter in eq. [I] so our 
results do not depend on the actual value of A w . 

The angular correlation function depends on UV lumi- 
nosity (see e.g. Giavalisco &; Dickinson 2001 Heinis et al. 
2007||Savoy et al.||2011| ), while there" is no evidence that it 
depends on f3. We checked that while there are some varia- 
tions in amplitude and slope of the correlation function with 
UV luminosity, they are not sufficient in this context to sig- 
nificantly change the correction for the flux determination. 
We consider only the best fit to the correlation function from 
the full sample hereafter. 

We show as triangles on Fig.[3]the ratio of the flux densi- 
ties corrected from clustering to the flux densities measured 
with PSF-fitting, as a function of Lfuv- The amplitude of 
the correction is blue around —20 per cent compared to the 
PSF-estimated flux at Lfuv < 10 10 L©, and becomes negli- 
gible for Lfuv > 10 10 L©. The amplitude of this correction 
is larger for the faint bins, where we observe a stronger de- 
parture from a PSF profile (see Appendix [51. 
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Figure 3. Amplitude of the correction applied to stacked flux 
densities as a function of Lfuv> a t 250, 350 and 500 jum. We 
show here the ratio of the corrected stacked flux density to the 
flux density measured by a simple PSF-fitting. The triangles show 
this ratio for the correction due to the clustering of the input cat- 
alogue, the squares are the correction applied because of the in- 
completeness of the input catalogue, and the circles show the ratio 
for the total correction including clustering and incompleteness. 



3.3 Summary: flux density measurements 

For a given stacking measure, we correct first for stack- 
ing bias by subtracting the stacking bias profile from the 
stack profile, adding errors of the profiles in quadrature. The 
stacking bias profile used is derived by stacking fake sources 
recovered by source extraction using the same Lfuv distri- 
bution as the galaxies in the stack under study. We then 
fit the resulting profile with eq. |l]), leaving S and a as free 
parameters, in order to get a flux density corrected from 
the effect of clustering. We show in Fig. [3] the amplitude 



of the correction (combining effects of incompleteness and 
clustering as circles) as a function of Lfuv- These correc- 
tions partly compensate for each other; the amplitude of the 
overall correction decreases with luminosity. The correction 
is maximal at Lfuv = 6 x 10 9 Lq and becomes negligible for 
luminosities larger than 2 x 10 10 Lq. The overall correction 
is larger at faint luminosities because the incompleteness is 
more important, and also because of the larger departure of 
the observed profiles from a pure PSF. Note that the errors 
on Fig. [3] are large for faint luminosities because the errors 
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on Spsf are large (around twice the errors on *Scorr ): the 
amplitude of the correction itself is well determined. 

The errors on the flux density S are obtained by boost- 
rap resampling, repeating the above procedure on 3000 ran- 
dom bootstrap samples, and determining the error from the 
standard deviation of the fluxes of these bootstrap samples. 

We further performed the following tests to assess the 
reliability of our results. First of all, we performed the full 
analysis using median stacking, and obtained very similar 



results (see Appendix CI ). We also tested how detected ob- 
jects impact our stacking results. We subtracted from the 
SPIRE images the sources detected at various threshold lev- 
els (from 3 to 10cr). We performed the stacking on these im- 
ages, and added to the flux measured by stacking the flux of 
the detected objects. The results we obtain with this method 



are slightly higher (on average 20 per cent, see Appendix C2 ) 



than the results presented here, while they agree at the la 
level. The results do not depend significantly on the thresh- 
old level we use to subtract detected objects. We discuss in 
Sect. |5.3| the impact of this difference on the cosmic star 
formation rate density derived from the IR luminosity func- 
tions we build from the stacking measurements. 



4 RESULTS 

We perform stacking at 250, 350, and 500 /im as a func- 
tion of UV luminosity and slope of the UV continuum, j5. 
Tables with stacking results are in Appendix [A] and we 
show the postage stamp images of the stacking in Appendix 
[B] We derive total IR luminosity, Lir, by integrating over 
the range 8 < A < 1000 /xm of the best fit to the |Dale fe| 



Helou (2002) templates, obtained with the SED-fitting code 
CIGALE pToll et al.|2009) . There are various ways to assign 
a redshift to a given stack population; we consider here as 
the redshift the mean photometric redshift of the galaxies 
involved in the stacking. The error on Lir is given by the 
standard deviation of the probability distribution function 
of the Lir values obtained with the models used during the 
fitting procedure (see Noll et aL|) 2009 ) . 

The Dale & Helou ( 2002 ) templates are calibrated as a 



function of a single parameter a. These models assume that 
the dust mass over interstellar radiation field ratio varies as 
a power law of the interstellar radiation field with index —a. 
We observe slight variations of this parameter as measured 
by the SED fitting. As a function of luminosity, we find 
that a decreases from ~ 2.6 for Lfuv < 1.5 x 1O 1O L to 
~ 1.7 for Lfuv > 1.5 x 10 10 Lq, which implies an increase 
of dust temperature from ~ 25 to 30 K. The significance of 
this increase is however only at a level of l.lcr. As a function 
of the slope of the UV continuum, /3, a is roughly constant 
at around 1.8, which implies a dust temperature of around 
28 K. 



4.1 Stacking as a function of UV luminosity 

We stack UV-selected galaxies in bins of UV luminosity 
Lfuv of size 0.1 dex. The results are presented in Fig. 




Figure 4. IR-to-UV luminosity ratio as a function of UV lu- 
minosity. Estimates with S/N > 3 at 250, 350 and 500 fim are 
shown as filled circles, and others as empty circles. The horizon- 
tally hatched region represents the locus of LIRGs, while the ver- 
tically hatched region is for ULIRGs. The dashed line shows the 
mean Lir/Lfuv value from estimates with S/N > 3 at 250, 350 
and 500 /im. The filled square shows the estimate from stacking 
of UV-selected galaxies at Herschel-PACS wavelengths at z ~ 2 
of |Reddy et al.| (2012) , slightly offset in Lfuv f° r clarity. The 
right axis shows the equivalent attenuation in the FUV band (in 
magnitudes), using eq. [2] The top panel shows the histogram of 
galaxies as a function of Lfuv- 



where we plot the ratio of IR to UV luminosities (which is a 
tracer of dust attenuation) as a function of the UV luminos- 
ity. The Lir/Lfuv ratio is found to be constant with UV 
luminosity over most of the range we probe, with a mean 
of 6.9 =b 1. This suggests that the dust attenuation does not 
depend heavily on UV luminosity in a UV-selected sample 
at z rsj 1.5 in this luminosity range. Assuming that the rela- 
tion between the attenuation in the GALEX FUV band and 
Lir/Lfuv is ( |Overzier et al.|2011| [Seibert et al.1|2005 ) 



A, 



2.5 log 



1 

L68 



+ 1 



(2) 



the average value of the IR to UV luminosity ratio for our 
sample corresponds to a value Afuv = 1.8 =b 0.1 mag. 

We also show in Fig. [4] the regions where Luminous In- 
frared Galaxies (LIRGs, 10 11 < Li R /L < 10 12 ) and Ultra 
Luminous Infrared Galaxies (ULIRGs, Lir > 10 12 Lq) lie. 
Our results show that the average IR luminosities of UV- 
selected galaxies in the UV luminosity range we explore at 
z ~ 1.5 are comparable to LIRGs, but not to ULIRGs. Our 
measures of the Lir/Lfuv ratio are in agreement with pre- 
vious measures obtained from objects within UV-selected 



4 We consider here as signal-to-noise ratio the ratio of the flux 
measured after applying the corrections described in sect.[3]to the 



error. Note that some stacking measurements can have observed 
signal-to-noise ratios (i.e. before applying any correction) lower 
than 3. 
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samples detected in both the UV and IR (Reddy et al. 2006b 



Buat et al. 2009 Reddy et al. 2010b. Our result is also in 



excellent agreement with the stacking study of Red dy et] 
al. (2012), who found that the average IR-to-UV luminos- 
ity ratio of a sample of 114 UV-selected galaxies at z ~ 2 
is 7.1 ± 1.1, for galaxies with L F uv ~ 3.1 x 10 10 Lq. Note 
also that Bu at et al.| ( 2009 ) observed that the fraction of 



galaxies with Lir/Lfuv > 5 is roughly constant with UV 
luminosity for Luv ^ 3 x 10 9 Lq, which is consistent with 
the trend we measure here. We can also compare our re- 
sult with the average attenuation derived by Cuc ciati et al.| 
(2012) from SED fitting of galaxies selected in the /-band 



(around 3000 A rest-frame at the mean redshift of our sam- 
ple) with u* g'r'i zJHK s photometry and spectroscopic red- 
shifts. They derive Apuv = 2.17 mag for the same redshift 
range (1.2 < z < 1.7) as our study, which is slightly larger 
than what we measure, though note that we select galaxies 
at shorter restframe wavelengths. 



4.2 Stacking as a function of UV slope, f3 

The slope of the UV continuum has been shown to corre- 



late with the dust attenuation within galaxies (e.g. |Calzetti,| 
Kinney, & Storchi-Bergmann 1994; Me urer, Heckman, fc| 



Calzetti 1999). The use of the /3 slope offers an estimate 
of the dust attenuation from the rest-frame UV, without 
requiring far-IR data or spectral lines diagnostics. Calibra- 
tions have been derived from spectro-photometric samples 



of starburst galaxies at low redshifts (e.g. Meurer, Heckman, 



&; Calzetti 1999 Overzier et al. 2011 ), and are routinely used 



to derive dust attenuation at various redshifts, in particular 
using slopes derived from rest-frame UV colours. 

We use the u* , V, intermediate (IA427, IA464, IA484, 
IA505, IA527, IA574, IA624, IA679, IA709, IA738, IA767, 
IA827) and narrow-band (NB711, NB816) filters to compute 
the slope f3. We adjust the photometry to a simple power-law 
SED, fx oc A^, over the restframe wavelength range 1200 < 
A < 3000 A. This means that there are at least 9 bands 
available for the measure of /3, and at the mean redshift of 
this sample, z = 1.43, there are 12 bands available. 

Figure [5] shows the distribution of (3 as a function 
of UV luminosity. The mean UV slope for our sample is 
(f3) — — 1.2 ± 0.6. We find that the average slope of the 
UV continuum is mostly independent of the UV luminosity, 
while the dispersion slightly decreases with increasing Lfuv . 
We compare these measures with the results of Bouwens 
et al.| ( |2009t obtained from [/-band dropouts from WFPC2 
(F300W) data at z ~ 2.5. While the two distributions over- 
lap at the la level, at a given UV luminosity the UV-selected 
galaxies are redder at lower redshifts. This is in agreement 
with a more global trend observed from 2 ^ 7 to 2; ~ 2.5 
(Bouwens et al. 2009). Note also that at higher redshifts, 



fainter galaxies tend to be bluer (with a significance of 5cr, 
Bou wens et al.|2009 ), while we do not observe such a trend at 
lower redshift. We note that our sample contains less than 1 
per cent of "quiescent" galaxies, according to the criterion of 
Ilbert et al. (2010) based on the restframe NUV — R colour. 



Moreover the lack of dependence of f3 with UV luminosity 
remains if we split galaxies in restframe NUV — R colour or 
specific star formation rate. 

We show on Fig. [6] the IR-to-UV luminosity ratio as a 



This work • 
oirwens et al. (2009) □ 
z~2.5 



cn -1 




10.5 
iog(L FUV [LJ) 



11.5 



Figure 5. UV continuum slope (3 as a function of UV luminosity. 
The large filled circles show the mean and dispersion in 15 bins 
of luminosity for Lpuv < 5 x 10 10 Lq. Open squares show the 
mean (3 measured for {/-band dropouts from WFPC2 (F300W) 
data at z ~ 2.5 (Bouwens et al. 2009]). 
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Figure 6. Lir/^fuv vs - P from stacking. Estimates with S/N > 
3 at 250, 350 and 500 jum are shown as filled circles, and others 
as empty circles. Lines show various £ir/£fuv — P relations: 
Calzetti et aT] ( |2000| dotte d); |Overzier et al.] (|2011| dash ed) ; |Hao| 
et al.| ( |201l7 "dot-dashed); |Munoz-Mateos et al.| ( |2009| dot-dot- 
dashed); and the relation expected for the SMC extinction curve 
(from [Overzier et al.|201l] long dashed). The solid line shows the 
best fit to our measurements (considering only estimates with 
S/N > 3), and the grey area the range of relations implied by the 
lcr errors on the parameters. The top panel shows the number of 
galaxies in each bin of (3. 
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function of /3. Note that the stacking results do not cover the 
same range in /3 as in Fig. [5] the signal is small for f3 < —1.5. 

We observe a good correlation between the UV slope 
and dust attenuation. We compare our results to several 

calibrations of the (Lir/Lfuv) P relation derived for 

local starburst galaxies ( | Calzetti et al. 2000; Overzier et al 
|2011 [ their 'total relation'). The difference between these two 



local relations comes from the fact that Overzier et al. (|2011 ) 



remeasured the UV photometry for the M eurer, Heckman,| 
& Calzetti (1999) sample, hence possibly including fewer 



starburst regions. 

Our results fall between the local starburst and the lo- 
cal calibration for normal star-forming galaxies from |Munoz-| 



Mateos et al. (2009) and the relation expected for the Small 



Magellanic Cloud (SMC) extinction curve. Our results are, 
however, consistent with the relation obtained by |Hao et al.| 
(2011 ) from another set of local normal star-forming galax- 
ies, based on SINGS observations ( Kennicutt et al.|2003 ), as 



well as data from Moustakas & Kennicutt 



the spread between the calibrations from 
al 



(2009 ) and Hao et al. (2011) comes from differences in 



!2006). 



Note that 



Muhoz- Mateos et] 



galaxy selections. 

Following previous work we fit our measurements as- 
suming 



2.5 log 



1 

L68 



+ 1 



: a + bfi. (3) 



We consider only measurements with S/N > 3 in the 
fit; including other measurements does not affect the results. 
We find a = 3.4 ± 0.1 and b = 1.6 ± 0.1. Our value for the 
slope b = dAFuv/d/3 of this relation is lower than the val- 
ues derived from commonly used local starburst attenuation 
laws (e.g.|Meurer, Heckman, fe Calzetti| ( |1999| ) find b = 1.99; 



Calzetti et al. ( 2000 ) find b = 2.31), but larger than the value 
derived at intermediate redshifts (1 < z < 2) by|Buat et al. 



( 20111 , namely b = 1.46±0.21. The implied UV slope for the 
dust-free case is fto — — 2.12±0.18; this is in agreement with 



what is expected from stellar population models (Leitherer 
&; Heckman|19 95 ) , and favours a continuous star- formation 



mode. 



5 UV AND IR LUMINOSITY FUNCTIONS 

In this section, we present our determination of the UV lu- 
minosity function of our UV sample. Using the stacking re- 
sults presented above, we can also recover the total far-IR 
luminosity function of this sample. This procedure enables 
us, for instance, to discuss the dust-corrected contribution 
to the star- formation rate density of UV-selected galaxies. 

Reddy et al.| ( |2008| performed a similar study on a sam- 
ple of LBGs at 1.9 < z < 3.4. They derived the E(B - V) 



distribution of their sample by maximising the likelihood 
of observing their data for a given luminosity, redshift and 
reddening distribution. Assuming the |Meurer, Heckman, fe| 
Calzetti (1999} attenuation law, they determined the IR lu- 
minosity function of UV-selected galaxies using a Monte 
Carlo method and found that at z ~ 2 it is in agreement with 
the IR luminosity function of 8 /im rest-frame selected galax- 
ies with luminosities 10 10 — 1 12 Lq range . This contrasts 



with the lower redshift result of Buat et al. (2009), who di- 
rectly measured the IR luminosity function of UV-selected 



galaxies at z ~ 1 and noticed that it underestimates the IR 
luminosity function of 12 /im rest-frame selected galaxies for 
Lir>2x 10 n L o . 

We derive luminosity functions using the Vmax method 
(Schmidt 1968). In practice, we derive for each galaxy of 
our sample the minimum (z m in) and maximum (z max ) red- 
shifts where it can be included in the sample given its red- 
shift and luminosity: z m i n = max(1.2, zuv,min), and z ma x = 
min(1.7, zuv,max). ^uv,min and zuv,max are the minimum 
and maximum redshifts implied by the magnitude limits we 
used to build our UV-selected sample. 

The maximum volume within which this galaxy can be 
observed is then given by: 



V(z) 



= V(z n 



V(Z U 



(4) 



A 3/ \ 
^ r cO), 



where A is the solid angle covered by the observations, and 
r c (z) is the comoving distance. We corrected for incomplete- 
ness as a function of luminosity using the simulations de- 
scribed Sect. [2] We define our completeness limit as the lumi- 
nosity where the incompleteness is equal to 20 per cent; this 
corresponds to Lfuv = 5 x 10 9 Lq and Lir = 5.6 x 10 10 Lq. 
We also include the error on the completeness correction 
in the luminosity function errors. For each luminosity func- 
tion estimate, we take into account the error on photometric 
redshifts by constructing 50 mock catalogues with new red- 
shifts within the probability distribution functions derived 



by Ilbert et al. (2009|. Note that this procedure yields an 
estimate of the errors added by the use of photometric red- 
shifts, but not of the bias they introduce. 

5.1 UV luminosity function 

We checked that the K-corrections are minimal and do not 
have an impact on the UV luminosity function. We note that 
the faint-end of the UV LF is quite sensitive to the method 
used for the photometry. The UV LF derived using the u*- 
band fluxes from the catalog of Capak et al.| p007) is signifi- 
cantly steeper (3a significance level, a — —1.29 ±0.03) at the 
faint end than the measurement obtained using the photom- 
etry we use here. [Capak et aL ( 2007| performed the source 
extraction in a combined i + and i* image, ran PSF match- 
ing to the image with worst seeing (K s ), and finally mea- 
sured aperture photometry and aperture corrections. We at- 
tribute this difference between the LFs to an overestimation 
of the u* flux because of a combination of the PSF match- 
ing and the aperture corrections. We show our UV LF in 
Fig. [7] (blue circles). The error bars on the UV LF are the 
combination of the Poisson error, the error on the complete- 
ness correction, and the standard deviation of the estimates 
from the mock catalogues used to compute the LF. We fit 
the UV luminosity function with a Schechter form; we find 
logO*[Mp C - 3 ]) = -2.24 ±0.02, log(L*[L ]) = 10.13 ±0.02 
(equivalent to M* = -19.44 ± 0.04), and a = -1.13 ± 0.04. 
We note that the luminosity function we derive is quite 
flat with respect to a number of previous studies devoted 
to estimating the UV luminosity functions at similar red- 
shifts ( Arnouts et al.|[2005| |Oesch et al.||2010 ), who report 
— 1.6 < a < —1.4. This difference could be due to the fact 
that these estimates were obtained from selections using 
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Figure 7. UV and IR luminosity functions. The blue circles rep- 
resent the UV luminosity function of our sample. The red circles 
show the IR luminosity function of this sample derived from the 
stacking measurement of the IR to UV luminosity ratio: filled cir- 
cles using a constant dispersion of log(LiR/Lpuv) (method '<x c '); 
open thick circles the dispersion which reproduces the observed 
log(LiR/Lpuv) values for detected objects ('cr^'); and open thin 
circles the dispersion derived from the dispersion in (3 slope ('cr«'). 
The dashed line shows the lower luminosity limit that we consider 
to adjust the IR luminosity function. Open squares show the IR 
luminosity function of an IR selected sample at 1.3 < z < 1.8 
from Ma gnelli et al.| (|2011|); op en triangles the IR luminosity 
function from [Rodighiero et al. (2010) at 1.2 < z < 1.7; and 
open hexagons the IR luminosity function from Herschel/ PACS 
data at 1.2 < z < 1.7 ( |Gruppioni et al.|2010] >. 



shorter restframe wavelengths than ours. However, our es- 
timate is similar to a recent and independent derivation of 
the UV luminosity function by Cucciati et al. (2012), who 
found a = —1.09 =b 0.23 in the same redshift range, using 
spectroscopic data. 



5.2 Recovering the IR luminosity function 

From the stacking results presented above, we recover the 
IR luminosity function of our UV-selected sample. To do so, 
we assign an IR luminosity to each of the galaxies of our 
sample using the following method, which is similar to that 
used by |Reddy et al.| ( |2008 ). 

We assume that at a given UV luminosity, the distribu- 
tion of log(LiR/LFuv) is Gaussian. We use as the mean of 
this distribution the values of Lir/Lfuv derived from the 
stacking analysis. In the luminosity range where we do not 
have reliable stacking measurements (Lfuv > 5 x 10 10 Lq), 
we assume that the Lir/Lfuv ratio is constant, and equal to 
6.8, the average value of the ratio over 10 10 < Lfuv/L© < 
5 x 10 10 . We also consider as mean of the Lir/Lfuv dis- 
tribution the results obtained from an alternative stacking 
method, excluding sources detected at 10cr from the SPIRE 



images (see Appendix C2). 

For the standard deviation of the Gaussian distribution, 
we use three different methods. Firstly, (i) we assume that 
the dispersion of the Gaussian distribution is constant with 
UV luminosity and equal to 0.35. This value has been de- 



rived at low redshifts from UV-selected samples by Buat et 
al. (2009); we call this method '<r c '. Alternatively, (zz), we use 
as reference the values of the Lir/Lfuv ratios of the objects 
detected at Herschel-SPIRE wavelengths (see Sect. [2]); these 
objects represent the upper tail of the distribution of the 
Lir/Lfuv ratio. We determined the mean of this distribu- 
tion from stacking as a function of Lfuv (Sect. 4.1 ). We as- 



sume once again that the distribution of the Lir/Lfuv ratio 
is Gaussian. Then for each Lfuv stacking bin, we adjust the 
standard deviation of the Gaussian such that we reproduce 
the observed distribution, given the mean and the few detec- 
tions in the upper tail. The resulting dispersion is a function 
of Lfuv, and decreases from 0.5 at Lfuv = 3 x 1O 9 L0 to 0.2 
at Lfuv = 10 1:l Lq (method 'era')- Finally, (Hi) we use the 
observed dispersion of the slope of the UV continuum, /3, as a 
function of UV luminosity, and translate it into a dispersion 
in log(LiR/LFuv), using the relation we derived in Sect. 4.2 
Apuv = 3.4+ 1.6/3. This also yields a function of Lfuv; the 
resulting dispersion is slightly higher than the one obtained 
with scenario cr d , decreasing from 0.6 at Lfuv = 3 x 10 9 Lq 



lO n L (method 'ap') 



to 0.3 at Lfu 

For a given galaxy in our sample, we then randomly as- 
sign a value of log (Lir/Lfuv), following the relevant distri- 
bution, whose mean and standard deviation are determined 
by the UV luminosity of the galaxy. We can then derive 
the IR luminosity for each galaxy in our sample, and, using 
the Vmax values determined according to the UV selection, 
compute the IR LF of the sample. 

At z = 0, all UV-selected galaxies are detected in the 
IR, and the distribution of log (Lir/Lfuv) is well described 
by a Gaussian ( |Buat et al.|20 09 ) . Note that the actual distri- 
bution of log(LiR/LFuv) at z ~ 1.5 is not known. Assuming 
a Gaussian distribution enables us to compare with previ- 
ous studies similar to ours (e.g. |Reddy et al.|2008| ). However, 
depending on the value of the dispersion of this Gaussian, 
a significant number of galaxies with low Lfuv can be as- 
signed high Lir. 

We generate 100 realisations of this IR LF, and show 
on Fig. [7] the mean and errors of these iterations (red cir- 
cles). We compare our result with the IR LF derived at 
1.3 < z < 1.8 by Magnelli et al.| (|2011| from a sample of 



galaxies selected at 24 /xm and also using stacking at 70 /xm 
(open squares), the IR LF of Rodig hiero et al.| ( |2010| (open 
triangles), from a sample of galaxies at 1.2 < z < 1.7, also 
selected at 24 /xm with mid-IR data, and finally the Her- 
schel /PACS- derived IR LF from |Gruppioni et al.| ( |2010| ) at 
1.2 < z < 1.7 (open hexagons J] 

Our results show that correcting a UV-selected sample 
for dust enables us to recover the IR LF at the faint lumi- 
nosities reached by the IR selections (Lir < 3 x 10 11 Lq); 
our estimates are in agreement with the results of Magnelli 
|et al.| ( |2ull| at these luminosities. 

At IR luminosities brighter than Lir = 3xlO n L , the 
results depend on the assumptions about the shape of the 
distribution of the Lir/Lfuv ratio. A higher dispersion in 
this ratio yields a higher amplitude of the LF at the bright 
end. We note however that for the scenario with the highest 
dispersion, erg, it is not clear that the dispersion in ft is com- 



5 We show in an Appendix (Fig. |C3| the IR LF obtained using 
the alternate stacking measurements. 
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pletely related to the dispersion in dust attenuation. First 
of all, we do not observe a stacking detection at S/N > 3 in 
all Herschel-SPIRE bands for /3 < -1.2. Note also that the 
dispersion in Lir/Lfuv that we infer from the dispersion in 
f3 is larger than the dispersion we derive from the objects 
directly detected by SPIRE. In any case, none of the scenar- 
ios we explore here for the dispersion of this ratio is able to 
reproduce accurately the bright end of the IR selected LFs. 
This is implied directly by the stacking results (see Fig. |4j), 
which show that the IR luminosities of our sample galaxies 
are consistent with LIRGs, but not ULIRGs. 

Hereafter, we consider as our best estimate the IR LF 
determined with the scenario <jd described above. 



5.3 Implications for Cosmic star-formation 
density estimation 

We study here the implications of our results for the estima- 
tion of the cosmic star-formation rate density (SFRD) from 
UV-selected samples. 

We compute the UV luminosity density using p*,uv = 
0*L*r(a + 2), and then convert it to a star formation rate 



density using the relation from Kennicutt (1998): 



SFRIMqyi- 1 } = 1.4 x 10" 28 L,[erg s" 1 Hz" 1 ], (5) 

which assumes a Salpeter ( |1955| ) Initial Mass Function 
(IMF). We find a star-formation density p*,uv — 0.022 ± 
0.001 M yr _1 Mpc" 3 ; the error on p*,uv is derived from the 
allowed range of UV LFs within la of the best fit. 

IR luminosity function is commonly described by a dou- 
ble power law (DPL, |Sanders et aL|2003) : 

<f)(L) = 0* f - — ) L < Lknee 



0(L) 



L 

-^knee 



L> L k 



(6) 



Alternatively one can use a double exponential (DE, 
Saunders et~aL|l990| ): 



0(L) = 



exp 



1 1 2 

-~2 log 



1 + 



We use both parameterisations, in order to compare 



with the results of Magnelli et al. (2011), who use the DPL 



form, and Rodighiero et al. (2010 ) who use the DE form. We 
allow to be free parameters the normalisations, 0knee and 
(/>*, the characteristic luminosities, Lknee and L*, and the 
bright-end slopes, a 2 and a. We keep the faint-end slopes 



fixed, using the same assumptions as Magnelli et al. (2011) 



(ai = -0.6) and |Rodighiero et al.| (2010| ) (a 1.2). The 
best-fit parameters are given in Table]TJfor our baseline IR 
LF. 

We then integrate the IR LF within two luminosity 
ranges: 10 7 < Lir/Lq < 10 15 , which implies an extrap- 
olation of the measured luminosity function to faint lumi- 
nosities, and 10 11 < Lir/Lq < 10 13 , which is the range 
where observations are available. We finally convert the IR 
luminosity density to a SFRD using 



SFRIMqjt- 1 } = 4.5 x 10" 44 Li R [erg s" 1 ] 



(8) 
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Figure 8. Cosmic star- format ion rate history. Our results are 
represented as filled circles: UV selection (blue); and IR LF of 
the UV selection (red). We compare these results to the estimates 
based on UV-selected samples (blue-like colours): Gab asch et al.| 
(open hexagons, |200 4); Schim movich et ah] (open circles, 2005| >; 
Wyd er et al.| (open square, 2005); Sawicki & Thompson (open 
stars, |2006| ; and |Cucciati et al.| (open triangles, |2012| ). We also 
show estimates based on IR selections (red-like colours) : |Le Floc'"h| 
et al.|(|2005| open hexagons) ; |Caputi et al.| ( |2007| open hexa gons) ; 
Rodighiero et al. (2010| open triangles); and |Magnelli et aT] ( |2011| 
open squares). 



(Kennicutt 1998), which also assumes a Salpeter IMF. The 
resulting star- formation rate densities are listed in Table [I] 
We consider first our most secure estimate of the SFRD, 
from the luminosity range where measurements are available 
(10 11 < Lir/Lq < 10 13 ). The SFRD derived from our base- 
line IR LF of the UV-selected galaxies is about 89 =b 32 per 
cent of the density derived from the result of Magn elTT et| 
SDJ20TT), and 52 ± 25 per cent of the density derived from 
Rodighiero et al. (2010) from IR-selected samples. If we con- 



sider the IR LF built from the alternative stacking results, 
these percentages are 112 zb 41 and 65 =b 37 per cent. 

We show our results for the cosmic star-formation rate 
extrapolating the luminosity function to low luminosity in 
Fig. [8] along with other UV-selected and IR-selected mea- 
surements (references on the figure), all integrated over the 
same luminosity ranges, and converted to the same IMF. 
The result from our UV selection is in good agreement with 
previous determinations of the star-formation rate density 
from UV-selected samples. 

The cosmic SFR we derive from the IR LF of the UV 
selection is in broad agreement with the estimates from IR- 
selected samples at similar redshifts. Nevertheless, given the 
differences in the estimates of the SFRD from IR-selected 
samples, the percentage of the IR SFRD recovered by our 
method varies from 100 ± 33 per cent if we co nsider |Mag- 



nelli et al. (2011 ) to 61 ± 27 per cent considering Rodighiero 
et al. (2010). If we consider the IR LF built from the al- 
ternative stacking method, these percentages are 131 ± 40 
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Table 1. LF fit parameters and derived Cosmic Star- Formation rate density 
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IR LF from UV, DPL fit IR LF from UV, DE fit 



logOknee) 




-2.70 ±0.02 


log(</>*) 


-2.31 ±0.04 


log(L knee ) 




11.49 ±0.02 


log(L*) 


11.10 ±0.09 


CX-2 




-2.10 ±0.08 


cr 


0.41 ±0.03 


ai 




-0.6 


a 


1.2 


P*,IR (10 11 


< L : r/L < 10 13 ) 


0.08±0.01 




0.08±0.02 


P*,IR (10 7 


< L : r/L < 10 15 ) 


0.16±0.01 




0.13±0.03 



DPL stands for Double Power Law, and DE for Double Exponential. ol\ 
and a are kept fixed during the fitting procedure. 0knee an d </>* are given in 
Mpc -3 , Lk nee and L* in Lq and p*,ir in MQyr _1 Mpc -3 . 



and 76 ± 38 per cent. This suggests that, after correction 
for dust attenuation, a UV-selected sample at z ~ 1.5 down 
to Lfuv = 3 x 10 9 Lq can recover the total SFRD esti- 
mated from the IR, the remaining uncertainty being on the 
discrepancy between the IR-selected LFs. 

Our results also imply that the dust-corrected estimate 
of the SFR density is roughly 6 times higher than the deter- 
mination from direct UV observations. 

We can compare our SFRD estimate with what 
would be derived using an average correction derived from 
f3. Specifically, we derive the average attenuation factor 
(SFRir + SFRuv)/SFRuv ^ 0.64Lir/L FU v + 1 from the 
distribution of following the same method as Bouw ens et| 
al. (2012). We use our best fit for the relation between Afuv 



and f3 (see Sect. 4.2). We can then determine the UV SFRD 
obtained from the UV LF integrated over the full data range 
(i.e. down to 0.2 x L*,fuv), and correct it for dust attenua- 
tion with this average factor. We compare this value to the 
SFRD derived from the integration of the IR LF of the UV 
selection over the full IR luminosity range we probe. The 
UV SFRD corrected for dust attenuation using f3 is around 
25 per cent lower than the SFRD derived from the IR LF 
of the UV selection. This suggests that, at z = 1.5 at least, 
using an average dust correction factor derived from /3 can 
lead to a significant underestimation of the SFRD. 



6 DISCUSSION 

6.1 Comparison with previous studies 

There have been a number of studies exploring similar topics 
at lower and higher redshifts than our selection. We compare 
these with our results for the relation between the slope of 
the UV continuum and the dust attenuation, as well as the 
measure of the IR LF of a UV-selected sample. 



6.1.1 Relation f3 slope - Lir/Lfuv 

We find that there is a correlation between the slope of the 
UV continuum and the ratio Lir/Lfuv- However, the re- 
lation we observe at z — 1.5 is different from the relation 
that is derived from calibrations performed at low redshifts 



on starburst galaxies (Meurer, Heckman, &; Calzetti 1999 



Calzetti et al. 2000). This starburst relation is commonly 



used over a wide redshift range to correct UV luminosi- 
ties for dust attenuation when far-IR measurements are not 



available (e.g. Adelberger h Steidel 2000 Schiminovich et 
aL|2005| [Bouwens et al.||2009| ) 



It is claimed that the Meurer, Heckman, fc Cafzetti| 



( 1999 ) and Calzetti et al. ( 2000 ) relations are valid at various 
redshifts: Reddy et al. ( |2006a| show that typical z ~ 2 UV- 
selected galaxies detected at 24 /xm do follow the Meurer et 
al. relation. [Magdis et al.| ( |2010a| ) and |Magdis et aI.|p0T0b ) 
also find that the UV star-formation rates of Lyman Break 
Galaxies at z ~ 3 corrected for attenuation with the Me urer, | 



Heckman, & Calzetti ( 1999|) relation are in agreement with 
far-IR and radio estimates. |Reddy et al. ( |2006a| find, how- 
ever, that the dust attenuation is overestimated for galaxies 
with younger stellar populations; they also show that galax- 
ies with lower SFRs tend to lie under the |Meurer, H eck- 
man, fc Calzettij ( |1999| ) relation. |Burgarella et al.| ( |2Q07| 
notice that UV star-formation rates corrected for dust at- 
tenuation with the Meurer et al. relation are in agreement 
with IR-derived ones, but that the dispersion is large. In a 
study similar to ours, Reddy et al. ( 2012| perform stacking 
of z rsj 2 Lyman Break Galaxies at 24, 100 and 160 /xm, as 
well as 1.4 GHz. They find that the average Lir/Lfuv ratio 
of these L* Lyman Break Galaxies is consistent with the 
local starburst relation. 

While our result hence might seem at odds with previ- 
ous work, a number of studies have shown that attempting 
to correct for dust attenuation using the slope of the UV 
continuum is quite complex, and requires some caution. 

At low redshift, iBelll 120021), iBoissier et al.l (120071), 



Munoz-Mateos et al. (2009), and Seibert et al. (20051 showed 



that the f3 slope and the dust attenuation of normal star- 
forming galaxies do not follow the same relationship as 
for local starbursts; in this case the IMeurer, Heckman, fc| 



Calzettij ( |1999| ) relation overpredicts the dust attenuation, 
which is in agreement with our results. The validity of the 
starburst relation depends also on the sample selection: for 
local galaxies dust attenuation might be overestimated for 
a UV-selected sample, while underestimated for IR-selected 
samples ( |Buat et al.|2005||Seibert et al.|2005| ). On the other 



hand, at higher redshifts (0.66 < z < 2.6), Murphy et 



al. (2011) observed that the |Meurer, Heckman, &; Calzetti 



(1999) relation overpredicts the dust attenuation for 24 fim 



selected galaxies. Note also that Buat et al. ( 2010| found 
that at z < 0.3, galaxies selected at 250 /xm have dust at- 
tenuations between the local starbursts and the normal star 
forming relation from Boi ssier et al.| ( |2007| ) . The mean rela- 
tion that we derive here is actually in excellent agreement 
with their results. We note also that our sample is domi- 
nated by galaxies of IR luminosities similar to LIRGs. At 
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z ~ 1.5, these galaxies belong mostly to the 'main sequence' 
of star forming galaxies and are not in the starburst mode 
(Sar gent et al.|2012 ), which can partly explain why the rela- 
tion between f3 and the dust attenuation we obtain is similar 
to normal star forming galaxies at low redshifts. 

The variety of results presented above is linked to the 
fact that the UV and IR emissions can come from differ- 
ent regions in the galaxies, and also that attempting to use 
the slope of the UV continuum to correct for dust attenu- 
ation requires assumptions about the underlying extinction 
law and star-formation history. Using the Meurer, Heckman, 
fe Calzetti| ( |1999| ) or the |Calzetti et al.| ( |2QQQ| ) relations is 
consistent with using quite shallow extinction laws, such as 



persion (Lir /Lfuv = 4.7, cr(log(LiR/L F uv)) = 0.53, |Reddy| 
et al.|2006a| ). This average IR-to-UV ratio is lower than the 



the Calzetti, Kinney, & Storchi-Bergmann ( 1994 ) or Calzetti 



(1997) ones, which is equivalent to adopting a clumpy fore- 



ground distribution of dust. Our results (in terms of the 
slope of the Afuv — ft relation) suggest a steeper extinction 
law, which is expected for dust geometries similar to a fore- 
ground screen. Note however that using the SMC extinction 
law yields a relation lower than our results, and hence is too 
steep. On the other hand, the value of the dust-free slope 
we derive suggests that the star formation mode is more 
continuous rather than starburst. 

In summary, our result for the f3 slope - Lir/Lfuv rela- 
tion at = 1.5 is similar to those obtained from local normal 
star forming galaxies. This suggests that our sample galaxies 
have a steeper extinction law than tho se from|Calzetti, Kin-"] 
ney, & Storchi-Bergmann| ( |1994| ) and |Calzetti| ( |1997| , but 
less steep than the SMC one. This also shows that using 
a /3 slope - Lir/Lfuv relation calibrated on local starburts 
galaxies can induce an overestimation of the Lir / Lfuv from 
the /3 slope by around a factor of 2. 

6.1.2 IR LF of UV- selected sample 

We find that the IR LF of a UV-selected sample at z ~ 1.5 
down to Lfuv = 3 x 10 9 recovers the faint-end of the IR 
LF of far-IR selected samples, but might underestimate the 
bright-end, if we consider the latest Herschel results, and 
our most conservative stacking measures. 

These results are in agreement with those of |Buat et] 



al. (2009) (z ~ 1) who found that the bolometric (L = 



Luv + Lir) LF of a UV selection directly measured from 
UV and IR data underestimates the IR LF from I R selected 
sample for L > 2 x 10 11 Lq. At higher redshifts, Reddy et 
|al.| ([2008) performed a similar study, and found, on the con- 
trary, that the reconstructed IR LF of UV-selected samples 
at z = 2 and z — 3 is similar to the IR LF of IR-selected 
samples. Note that we use a method which is quite similar 



to that of|Reddy et al.| ( |2QQ8| ) to build our IR LF. In detail, 
Red dy et al.| ( |2008| ) reconstruct two IR LFs: one from the 
distribution of E(B — V) derived from a maximum likeli- 
hood analysis; and the other one using previously observed 
Lir/Lfuv ratios and dispersion. These two LFs are consis- 
tent with each other. We note that their method of using the 
distribution of E(B — V) assumes the |Calzetti et a l. (2000) 
relation, which does not seem to apply to our sample (Sect. 
6.1.1). In other words, since using the Calzetti et al. (2000) 



relation on our sample would yield larger dust-corrected lu- 
minosities, this can explain part of the discrepancy. 

The other method used by Reddy et al.| ( |2008| is based 
on previous determinations of the IR-to-UV ratio and its dis- 



values we derive here, while the dispersion is higher. In or- 
der to test the impact of our assumptions on the recovery 
of the bright-end of the IR LF, we used as an extreme case 
the values of Lir/Lfuv we obtained here, and the dispersion 
used by |Reddy et al.| (2008) to construct the IR LF from the 
UV selection. The results we thus obtain (see filled squares 
on fig. [9| are in agreement with the IR LF of Magnelli et 
[aL ( 120111, based on Spitzer data), but then slightly under- 
estimate the measurement s of |Rodighiero et al. (2010 also 
based on Spitzer data) and |Gruppioni et al. (2010 based on 
Herschel PACS data). While this result is in better agree- 
ment with the LFs obtained from IR-selected samples, we 
believe that a constant dispersion of the (Lir/Lfuv) with 
Lfuv is unlikely. Based on our sample, we observe that the 
dispersion around (Lir/Lfuv) decreases with Lfuv, either 
from the dispersion in /3 or from the dispersion which is re- 
quired to reproduce the (Lir/Lfuv) values of the few UV- 
selected galaxies detected at the SPIRE wavelengths. Note 
also that using this value of the dispersion would imply, for 
instance, around four times as many ULIRGs as we detect 
with SPIRE. We also attempted to adjust the dispersion 
required to match the IR LF of Rodighiero et al. (2010), as- 



suming an exponentialy decreasing function of Lfuv, as we 
observe in the data. The best fit in this case implies twice 
as many ULIRGs as we detect. 

6.1.3 Contribution to the CIB of UV-selected sources 

Our results for the IR LF of our UV-selected sample imply 
that a UV-selection misses a significant part of the IR galaxy 
population, at least within the UV luminosity range we are 
able to probe. To investigate this further, we measured the 
contribution to the Cosmic Infrared Background (CIB) at 
the SPIRE wavelengths from the UV-selected galaxies in 
our sample, and compared it to the values measured by 
Bethermin et al. (2012), who derived the contribution to 



the CIB from 24 /xm-selected sources with $24 > 80 jjl Jy. We 
measured by stacking the average flux density for galaxies 
in two redshift bins (1.2 < z < 1.4 and 1.4 < z < 1.6), and 
converted it to a surface brightness. We find that the con- 
tribution to the CIB from our UV-selected galaxies is lower 
than that from 24/xm-selected sources. The contribution to 
the CIB from our UV-selected galaxies is around 50 per cent 
of that from 24 /xm-selected sources for 1.2 < z < 1.4 and 
around 30 per cent for 1.4 < z < 1.6. This result clearly 
shows that a UV selection is missing a part of the galaxy 
population probed by IR selections, and that the amount 
of energy which is emitted by this missing population is a 
significant fraction of the CIB. 

6.2 Recovering the bright-end of the IR LF 

Using our most conservative stacking measurements, our re- 
sults show that the IR LF we build for our UV-selected 
sample does not recover the bright-end of the measured IR 
LF from IR-selected samples, if we consider the latest Her- 
schel estimates. This would imply that a part of the IR 
galaxy population is missed by UV selection, at least in the 
redshift and luminosity ranges we consider here. We investi- 
gate the possibility that this missing population is actually 
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Figure 9. Recovering the bright-end of the IR LF. Filled circles 
show the IR LF of our UV selection, based on available data. 
The solid line show the total IR LF of the UV selection including 
the contribution of fainter galaxies (10 8 < Lfuv/L q < 10 9-5 ) 
adjusted to match the result of [Rodighiero ~et al.| ( |20To| . The 
dashed line show the contribution to the IR LF of the faint UV 
galaxies. The shaded grey area shows the range of LFs implied 
by the lcr errors on the parameters. We assume here that the 
dispersion of the ratio log(LiR/Lpuv) is a constant for galaxies 
with 10 8 < L FUV /L < 10 9 - 5 . The filled squares show the IR 
LF of the UV selection (no extrapolation to luminosities fainter 
than our limit) we obtain if we use our stacking measurements for 
log(LiK/LFuv) an d the dispersion on this ratio used by |Reddy| 
|et al.| ( [2008] ). 



fainter in UV luminosity than the limit reached by our ob- 
servations. Note that according to our results shown in Fig. 
[4] the galaxies in our UV-selected sample are mostly similar 
to LIRGs, hence we have to populate the ULIRGs regime of 
the luminosity function, which requires a large dispersion of 
the Lir/Lfuv ratio. 

To do so, we extrapolate the UV luminosity func- 
tion down to Lfuv = 10 8 Lq, using the best fit we ob- 



tain (see Sect. 5.1). This corresponds to a magnitude limit 
u* = 29.75 mag, 3.75 magnitudes deeper than the data we 
are using. In practice we create a mock catalogue which has 
the proper luminosity function. Then we compute the IR LF 
of this mock catalogue using the same method described in 
Sect. |5.2[ In this case, however, we do not have measure- 
ments for the mean IR-to-UV ratio or the dispersion of this 
ratio. We rather adjust these two quantities such that the 
sum of the IR LF of the mock catalogue and the IR LF 
of the UV selection from our baseline stacking method is 



consistent with the IR LF of Rodighiero et al. ( 2010| ) 



As we are limited by the relatively poor observational 
constraints on the IR LF, we make the following simplis- 
tic assumptions: we consider that the ratio Lir/Lfuv and 
the dispersion around this ratio are constant over the range 
10 8 < Lfuv/L© < 3 x 10 9 ; we then use two free parameters, 



(Lir/Lfuv) and a (log(LiR/LFuv)), to adjust the contribu- 
tion of UV faint galaxies to the bright-end of the IR LF. 

The results of the fit are shown Fig. [9] The contribu- 
tion to the IR LF of the faint UV galaxies is shown as the 
dashed line, and the total IR LF as the solid line. These 
results show that it is possible to obtain parameters that 
fit the IR LF measured from IR-selected samples. We find 
(Lir/Lfuv) = 8.7±i^ and a(log(Li R /L F uv)) = 0.73lg;gJ. 
The average IR to UV ratio required to match the bright- 
end of the IR LF is consistent with the value we mea- 
sure in our two faintest bins at Lfuv = (3 — 4) x 10 9 Lq 
(Lir/Lfuv = 14.5 =b 6, 11.5 =b 5), while the requested dis- 
persion of the ratio is larger than what we observe: at 
Lfuv = 3 x 10 9 Lq the dispersion required to match the 
detected objects is a(log(LiR/LFuv)) = 0.5. 

From this fit we can estimate the fraction of galaxies 
that is missed by a UV selection, compared to an IR selec- 
tion. We integrate the IR LF of the faint UV galaxies in the 
range 6.5 x 10 10 < Lir/L < 10 13 , where data are available, 
and compare this density to the one derived by integrating 
the total IR LF in the same range. We find that in terms of 
number density, 56 per cent of galaxies are missed. 

The relatively poor constraints on the IR LF do not 
enable us to draw firm conclusions on the scenarios of dust 
attenuation for faint UV galaxies. Note that the shape of the 
IR LF is quite sensitive to the assumptions on the dispersion 
a (log (Lir/Lfuv)) - Hence we expect to set better constraints 
on these scenarios with updated Herschel LFs. 



7 SUMMARY AND CONCLUSIONS 



We used Herschel-SPIRE ( jGriffin et aL]|2010| |Swinyard et 



al. 2010) imaging from the HerMES (Oliver et al. 2012) 



programme to study the IR properties of a sample of UV- 
selected galaxies at z ~ 1.5 in the COSMOS field. We built 



our sample from galaxies detected in CFHT u* band ( Capak 
et al. 2007) down to u* — 26 mag, with photometric redshifts 
plbert et aT||2009 ) 1.2 < z pYlot < 1.7. Only a few per cent 
of these galaxies are detected at the Herschel- SPIRE wave- 
lengths, so we use stacking in order to derive the average 
IR luminosities of the galaxies as a function of Lfuv and 
the slope of the UV continuum /3. We detail the techniques 
we used to correct the stacking measurements from stack- 
ing bias and clustering of the UV-selected galaxies, based on 
extensive simulations. We use these stacking measurements 
to derive the IR LF of the UV-selected sample, in order to 
infer the cosmic star formation rate density probed by a UV 
selection at z = 1.5. Our conclusions are as follows: 



1. UV-selected galaxies at z 



1.5 and 4 x 10 9 < 



^FUV /L < 5 x 10 10 have average total IR luminosities 
similar to LIRGs, but not to ULIRGs. 

2. The average Lir/Lfuv ratio is roughly constant, with 
Lfuv (4 x 10 9 < L F uv/L© < 5 x 10 10 ) and is equal to 
6.9 ±1. 

3. The average Lir/Lfuv ratio is correlated with the slope 
of the UV continuum, f3. This relation is below the re- 
lation derived from local starburst galaxies, but is in 
agreement with previous results obtained from local nor- 
mal star-forming galaxies. Our best fit to this relation is 
Afuv = 3.4 ± 0.1 + (1.6 ± 0.1)£. 
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4. We built the IR LF of the UV sample using our stacking 
measurements of the average Lir/Lfuv ratio, and assum- 
ing that the distribution of log(LiR/LFUv) is Gaussian. 
We used three different scenarios for the value of the dis- 
persion cr(log(LiR/LFuv)), which all yield the same result 
that the IR LF of the UV sample is in reasonable agree- 
ment at the faint-end (Lir ~ 10 11 Lq) with the IR LF 
from IR-selected samples at the same epoch, but might 
underestimate it at the bright-end (Lir >5 x 10 11 Lq). 

5. At z rsj 1.5 a UV rest-frame selection without dust at- 
tenuation correction probes roughly 10 per cent of the 
total (UV+IR) star-formation rate density. The cosmic 
star-formation rate density derived from the IR LF of 
the UV sample corresponds to 61-76 per cent or 100-133 
per cent of the star-formation rate density derived from 
IR-selected samples, depending on the IR LF taken as 
reference (Rodighi ero et al.|2010 or Magnelli et al.|20TT ). 

6. Assuming our most conservative measures and the lat- 
est Herschel estimates, the fraction of galaxies which are 
missed by a UV selection compared to an IR selection at 
z ~ 1.5 is around 50 per cent, in terms of number den- 
sity; this number is sensitive to the assumptions on the 
dispersion cr(log(Li R /L F uv)). 
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Table Al. Stacking results as a function of Lfuv 



log(L FUV [L©]) range (log(L FUV [£©])) (z) N gat S 2 50 [mjy] #350 [mjy] #500 [mJy] log(L m [L ]) 



9.44 


- 9.54 


9.51 


1 


.26 


872 


0.71 


± 


0.24 


0.84 


± 


0.21 


0.81 


± 


0.27 


10.67 


± 


0.20 


9.54 


- 9.64 


9.60 


1 


.34 


3075 


0.66 


± 


0.12 


0.72 


± 


0.16 


0.52 


± 


0.15 


10.66 




0.18 


9.64 


- 9.74 


9.70 


1 


.43 


5760 


0.46 


± 


0.08 


0.58 


± 


0.10 


0.44 


± 


0.10 


10.54 


± 


0.11 


9.74 


- 9.84 


9.79 


1 


.46 


6018 


0.43 


± 


0.10 


0.53 


± 


0.10 


0.36 


± 


0.09 


10.54 


± 


0.17 


9.84 


- 9.94 


9.89 


1 


.44 


5383 


0.63 


± 


0.09 


0.69 


± 


0.11 


0.54 


± 


0.10 


10.67 


± 


0.12 


9.94 - 


- 10.04 


9.99 


1 


.44 


4631 


0.65 


± 


0.10 


0.66 


± 


0.12 


0.48 


± 


0.10 


10.71 


± 


0.18 


10.04 


- 10.14 


10.09 


1 


.44 


3883 


0.88 


± 


0.11 


0.89 


± 


0.13 


0.59 


± 


0.11 


10.86 


± 


0.17 


10.14 


- 10.24 


10.19 


1 


.43 


3052 


1.35 


± 


0.13 


1.40 


± 


0.16 


1.02 


± 


0.13 


10.99 


± 


0.10 


10.24 


- 10.34 


10.29 


1 


.43 


2177 


1.54 


± 


0.17 


1.42 


± 


0.16 


0.90 


± 


0.14 


11.14 


± 


0.18 


10.34 


- 10.44 


10.39 


1 


.43 


1503 


1.87 


± 


0.21 


1.79 


± 


0.20 


1.10 


± 


0.19 


11.22 


± 


0.19 


10.44 


- 10.54 


10.49 


1 


.44 


873 


2.43 


± 


0.25 


2.38 


± 


0.29 


1.60 


± 


0.26 


11.29 


± 


0.15 


10.54 


- 10.64 


10.59 


1 


.44 


463 


3.45 


± 


0.40 


3.32 


± 


0.43 


2.10 


± 


0.34 


11.47 


± 


0.19 


10.64 


- 10.74 


10.69 


1 


.44 


211 


3.66 


± 


0.61 


3.39 


± 


0.61 


1.99 


± 


0.54 


11.56 


± 


0.24 


10.74 


- 10.84 


10.78 


1 


.43 


111 


3.56 


± 


0.77 


3.09 


± 


0.85 


2.42 


± 


0.85 


11.50 


± 


0.31 


10.84 


- 10.94 


10.89 


1 


.45 


35 


3.85 


± 


1.62 


2.46 


± 


1.92 


0.00 


± 


0.75 


11.18 


± 


0.96 


10.94 


- 11.04 


10.99 


1 


.43 


16 


9.50 


± 


1.58 


7.50 


± 


1.97 


4.69 


± 


1.92 


11.99 


± 


0.36 


11.04 


- 11.14 


11.07 


1 


.52 


5 


3.97 


± 


2.36 


3.89 


± 


3.51 


4.50 


± 


4.69 


11.59 


± 


0.58 


11.14 


- 11.24 


11.19 


1 


.44 


5 


22.96 ± 


: 8.25 


22.23 ± 7.21 


17.63 ± 5.04 


12.28 


± 


0.25 
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Table A2. Stacking results as a function of f3 



(3 range {(3) (log(L FU v[L©])) (z) N gat S 2 50 [mjy] #350 [mJy] #500 [mJy] log(Li R [L©]) 



-2.50 - 


-2.30 


-2.39 


9.92 


1 


.46 


418 


0.00 


± 


0.03 


0.08 


± 


0.17 


0.17 


± 


0.23 


8.46 ± 0.72 


-2.30 - 


-2.10 


-2.19 


9.98 


1 


.46 


921 


0.00 


± 


0.00 


0.00 


± 


0.00 


0.00 


± 


0.00 


0.00 ± 0.00 


-2.10 - 


-1.90 


-1.99 


9.98 


1 


.45 


1877 


0.00 


± 


0.00 


0.00 


± 


0.00 


0.00 


± 


0.00 


0.00 ± 0.00 


-1.90 - 


-1.70 


-1.79 


10.04 


1 


.45 


3448 


0.00 


± 


0.00 


0.00 


± 


0.00 


0.00 


± 


0.00 


0.00 ± 0.00 


-1.70 - 


-1.50 


-1.60 


10.05 


1 


.45 


4984 


0.00 


± 


0.00 


0.00 


± 


0.00 


0.03 


± 


0.04 


6.25 ± 0.66 


-1.50 - 


-1.30 


-1.40 


10.07 


1 


.44 


6081 


0.19 


± 


0.02 


0.38 


± 


0.09 


0.28 


± 


0.10 


10.22 ± 0.07 


-1.30 - 


-1.10 


-1.20 


10.07 


1 


.43 


5725 


0.44 


± 


0.09 


0.52 


± 


0.11 


0.38 


± 


0.11 


10.55 ± 0.17 


-1.10 - 


-0.90 


-1.00 


10.05 


1 


.42 


4705 


1.11 


± 


0.10 


1.09 


± 


0.11 


0.84 


± 


0.10 


10.89 ± 0.10 


-0.90 - 


-0.70 


-0.81 


10.02 


1 


.41 


3321 


1.76 


± 


0.13 


1.79 


± 


0.12 


1.18 


± 


0.10 


11.11 ± 0.11 


-0.70 - 


-0.50 


-0.61 


9.99 


1 


.41 


2109 


2.31 


± 


0.17 


2.22 


± 


0.17 


1.44 


± 


0.15 


11.24 ± 0.12 


-0.50 - 


-0.30 


-0.41 


9.99 


1 


.40 


1325 


3.40 


± 


0.22 


3.08 


± 


0.24 


2.09 


± 


0.20 


11.41 ± 0.12 


-0.30 - 


-0.10 


-0.20 


9.96 


1 


.39 


899 


4.59 


± 


0.30 


4.10 


± 


0.32 


2.72 


± 


0.27 


11.54 ± 0.13 


-0.10 - 


0.10 


-0.01 


9.94 


1 


.38 


604 


4.88 


± 


0.38 


4.48 


± 


0.40 


2.88 


± 


0.33 


11.56 ± 0.14 


0.10 - 


0.30 


0.19 


9.87 


1 


.37 


411 


6.38 


± 


0.49 


5.57 


± 


0.45 


3.42 


± 


0.43 


11.72 ± 0.16 


0.30 - 


0.50 


0.39 


9.87 


1 


.36 


266 


6.18 


± 


0.79 


5.71 


± 


0.71 


3.74 


± 


0.54 


11.63 ± 0.19 


0.50 - 


0.70 


0.59 


9.83 


1 


.36 


164 


6.26 


± 


0.59 


5.48 


± 


0.65 


3.17 


± 


0.59 


11.75 ± 0.20 


0.70 - 


0.90 


0.79 


9.77 


1 


.34 


107 


8.44 


± 


1.05 


8.68 


± 


0.99 


5.69 


± 


0.72 


11.71 ± 0.13 
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Figure CI. IR-to-UV luminosity ratio as a function of UV lu- 
minosity. Circles show estimates using the method we describe in 
this paper. Squares show estimates using median stacking, cor- 
recting for clustering only. The horizontally hatched region rep- 
resents the locus of LIRGs, while the vertically hatched region is 
for ULIRGs. The right axis shows the equivalent attenuation in 
the FUV band (in magnitudes), using eq. [2] The top panel shows 
the histogram of galaxies as a function of Lpuv- 



APPENDIX B: STACKING POSTAGE STAMP 
IMAGES 

We display here only the stacking results with signal-to-noise 
ratios larger than 3 at 250, 350 and 500 /xm. This signal- 
to-noise ratio is the ratio of the flux density over the er- 
ror on the flux density after the corrections described in 
Sect. [3] have been applied. We show here the raw stacking 
images, before applying any correction. All stack images are 
240 arcsec accross. The gray-scale shows the signal, using a 
asinh stretch ( Lupton et al.|[2004| . Each panel uses a dif- 
ferent scale, driven by the maximum flux. We also show for 
reference the radial profile along with the fit by the proper 
PSF in each band. 



APPENDIX C: COMPARISON WITH OTHER 
STACKING METHODS 



Comparison of mean stacking versus median 
stacking 



CI 



In Fig. |Cl| we show a comparison of the results obtained with 
our method with those obtained using median stacking. This 
figure is similar to Fig.^] We show in black circles the results 
obtained with our method, and in orange squares the results 
obtained with median stacking. For the latter, we perform 
median stacking, correcting only for clustering as described 



in Sect. 3.2 but not correcting for stacking bias (Sect. 3.1). 
The agreement between the results from mean and median 
stacking is excellent. 



8000 
6000 
' 4000 
2000 
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* Subtraction at 1 0o 
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Figure C2. IR-to-UV luminosity ratio as a function of UV lumi- 
nosity. Squares (respectively triangles) show the mean stacking 
results obtained while subtracting from the images the sources 
detected at 3a (respectively 10cr) in each band, correcting from 
clustering only. Circles show the results obtained with the method 
described in this paper. The horizontally hatched region repre- 
sents the locus of LIRGs, while the vertically hatched region is 
for ULIRGs. The right axis shows the equivalent attenuation in 
the FUV band (in magnitudes), using eq.|2] The top panel shows 
the histogram of galaxies as a function of Lpuv- 



Comparison of stacking with and without 
subtraction of detected sources 



C2 



We compare in Fig. |C2| the results of our method with an- 
other approach where we subtract detected sources from the 
images prior stacking. In practice we perform stacking after 
subtracting from the images the sources detected either at 
3cr or 10cr in each SPIRE band. We use mean stacking, cor- 
recting only for clustering as described in Sect. |3.2] but not 



correcting for stacking bias (Sect. 3.1). To obtain the final 
estimate, we add to the stacked flux measure the flux of 
the UV-selected sources detected with SPIRE which were 
subtracted. The resulting fluxes are on average 20 per cent 
higher than those obtained with our main method, even if 
they are in agreement at the la level, and do not depend 
significantly on the threshold used for subtracting sources. 
We show in Fig. |C3| the IR LF built using the method ad 



see Sect. 5.2) and these stacking measurements, and com- 



pare it to the results obtained with our baseline stacking 
results. This IR LF has a slightly larger amplitude at the 
bright end than the one obtained from our baseline results, 
and is hence in better agreement with the LF of Magn eTTi et| 
aD([20TT). 



REFERENCES 

Adelberger K. L., Steidel C. C, 2000, ApJ, 544, 218 
Arnouts S., et al., 2005, ApJ, 619, L43 



© 0000 RAS, MNRAS 000, 000-000 



18 S. Heinis et al 





Figure Bl. Stacking postage stamp images at 250 /im in bins of Lfuv- N g is the number of UV-selected galaxies in each bin. 
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Figure B2. Stacking postage stamp images at 350 /im in bins of of Lpuv- 
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Figure B3. Stacking postage stamp images at 500 jum in bins of of Lpuv- 
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1.2N g =5725 



0.6 N a =2109 




Figure B5. Stacking postage stamp images at 350 /im in bins of (3, the UV slope. 
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L[Ld 

Figure C3. IR luminosity functions. The filled circles show 
the IR luminosity function of our sample derived from the 
stacking measurement of the IR to UV luminosity ratio: the 
black circles using the dispersion which reproduces the observed 
log(LiR/Lpuv) values for detected objects ( L (Jd) an d our baseline 
stacking measurements; the red circles using the same scenario for 
the dispersion cr^, but using the stacking measurements obtained 
while subtracting from the images the sources detected at 10cr in 
each band. Open squares show the IR luminosity function of an 
IR selected sample at 1.3 < z < 1.8 from |Magnelli et al.| (|2011|>; 
open triangles the IR luminosity function from |Rodighiero et al.j 
(2010) at 1.2 < z < 1.7; and open hexagons the IR luminosity 
function from Herschel /PACS data at 1.2 < z < 1.7 (Gruppioni 
|et al.|2010| >. 
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